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We present lattice results on the scalar D s meson and comment on the 1)^,(2317) state recently 
discovered by BaBar and confirmed by CLEO, in view of a series of theoretical claims and counter 
claims. Lattice predictions in the static limit indicate larger masses than observed for a scalar quark 
model state. Finite c quark mass corrections seem to further enlarge this discrepancy, in support of 
a non quark-antiquark-state interpretation of experiment. 
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I. INTRODUCTION 

Recently the BaBar Collaboration announced the dis- 
covery of a Df positive parity meson at 2317±3 MeV and 
a width smaller than 10 MeV in the D+(1969)7r° chan- 
nel 0. No spin assignment has been made as yet but in 
view of the low mass J = appears likely. For simplicity 
we shall refer to it as the scalar or the + D s meson. This 
state was subsequently confirmed by the CLEO Collabo- 
ration 2] . CLEO also reports the observation of another 
state near 2460 MeV in the D*+(2112)7r° channel which 
is consistent with having J p = 1 + . These discoveries 
triggered a series of articles with different claims. In this 
note we discuss the scalar state in view of recent lattice 
results, after briefly summarising the different interpre- 
tations. 



II. QUARK MODEL OR NOT? 

If one treats the charm quark as a heavy spectator, the 
spin and angular momentum of the light antiquark can 
either couple to j = | (I = 0) or to j = | + and j — | + 
(I = 1). The interaction with the spectator spin will then 

result in a pseudoscalar- vector mass splitting for j = \ ■ 
in + and 1 + states for j = i + and in 1 /+ and 2 + states 

for J = § • The two 1 + states can undergo mixing. The 
pseudoscalar and vector D s states have been identified 
as L>+(1969) and D*+(2112), respectively. Then there 
is a £^(2536) state and a Djj(2573) which, with the 
likely spin assignment J = 2 in the latter case, form the 
j — | doublet. The j = \ + states can strongly decay into 
DK and D*K and are expected to be broad resonances. 
The new D+ 7 (2317) and the state at 2.46 GeV might 
constitute the missing doublet, where at least the former 
state, which lies almost 40 MeV below the DK threshold, 
is narrow. Cahn and Jackson Q interpret experiment in 
this way, in the context of a potential model. 

Barnes and collaborators [4| in contrast argue that 
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this state is most likely a DK molecule since its mass is 
160 MeV lighter than other potential model predictions 
which result in a mass around 2.48 GeV [j| for the scalar 
P-wave D s . Their argument is supported by the proxim- 
ity to the DK threshold and they interpret this system as 
a generalisation of an ao//o(980) KK molecule. A four- 
quark interpretation is also shared by Cheng and Hou [|| . 
Szczepaniak [7J argues in favour of a strong D s ir atomic 
contribution. 

Van Beveren and Rupp also liken this state with 
the ao//o(980) but interpret it as a quark model state. 
In their view the ao//o states are part of a low lying 
scalar quark-antiquark nonet, together with a tr(600) and 
a k(800). Consequently, they postulate additional scalar 
D mesons. According to them, in both the ao//o and 
the new D~^j systems, due to mixing with the KK or 
the DK continuum, respectively, the lowest scalar nonet 
is artificially lowered with respect to the quark model 
expectation. 

Bardeen et al. M discuss the heavy quark limit. They 
then follow Refs. |lOl ITlT ] and interpret the 0~ - + split- 
ting in terms of chiral symmetry. The symmetry break- 
ing scale corresponds in leading order to the constituent 
quark mass in the chiral limit ^3] an d has been esti- 
mated to be [13 AM w 338 MeV, a value that is very 
close to the experimental splitting of w 349 MeV. Chiral 
loops however will somewhat reduce the former expec- 
tation ITo|| . Colangelo et al. .12] share this picture and 
Godfrey investigates the decays that one would ex- 
pect in the case of a quark-antiquark interpretation. 

One should note that the vector-scalar splitting, which 
vanishes in the heavy quark limit, is as large as 143 MeV 
in the D s system, indeed an (D(A/m c ) correction to AM. 
In view of this, we would not expect the static approxima- 
tion to be quantitatively correct for D systems. We also 
remark that the + can be interpreted as a chiral part- 
ner of the 0~, independent of the quark model content, 
as long as isospin and strangeness agree. Unfortunately, 
most predicted decay rates in many of the above pictures 
seem to be more dictated by the mass and quantum num- 
bers of the state than by its quark content. However, in 
the case of an interpretation as a molecule or as part 
of an additional low lying scalar nonet (or triplet), an 
extra quark-model scalar state should still exist above 
the DK threshold. However, this might turn out to be 



2 



a rather broad resonance [f|. In contrast, in a straight 
D s interpretation there is no room for extra states other 
than the £>£(2460) and L»+ 1 (2536) between a £^(2573) 
(J = 21) and the newly discovered D+ J (2317) (J = 0?). 
The chiral heavy quark inter pre tation results in similar 
predictions for B systems 0, Hfl LUl j which in principle 
can be checked experimentally. 



III. THE STATIC LIMIT 

We will confront the new scalar state with lattice re- 
sults in the static limit in the quenched approximation 
and for rif = 2, before discussing finite charm quark mass 
corrections. 

In the static limit the j = \ and j = | doublets will 
be exactly mass degenerate. We wish to calculate the 

V and \ + 



masses. These can be extracted from the 
asymptotic large t decay of the two Euclidean correlation 
functions, 
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respectively. We made use of the relation M< = 75M75 
for the Wilson-Dirac operator M and 7475 = —7574. Uo,t 
denotes the Wilson-Schwinger line, connecting the point 
(x, t) with (x, 0). The spatial coordinate x is suppressed 
in U as well as in M and the colour trace is implicit. 
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FIG. 1: The 1 — i splitting in the static limit, as a func- 
tion of the light quark mass tx mp/my. Open symbols denote 
inter- and extrapolations to physical up/down and strange 
quark masses. Circles denote the nf — 2 case while squares 
denote the quenched approximation. The horizontal lines are 
the experimental values for the Z)j~j (2317) - Z)j~(1969) split- 
tings and the DK threshold. 

This static-light splitting has been calculated in the 
quenched approximation by Michael and Peisa |T3 | with 



Wilson action at (3 = 5.7 and (3 = 6.0. The results are 
in agreement with earlier references 0, 0, at addi- 
tional lattice spacings and no significant lattice spacing 
dependence has been observed: the fine structure split- 
tings, that strongly depend on short distance physics , 
vanish by definition. The quenched results, extrapolated 
to up/down and strange quark masses are depicted in 
Figure n (open squares) where r^ 1 w 400 MeV. The re- 
sults are AM S = 384(50) MeV and AM U = 299(114) 

MeV, respectively. The splittings of the i states with 

respect to the | + states are 434(5) MeV and 323(131) 
MeV for strange and up/down quarks, roughly 50 MeV 
larger. In contrast, the -0^(2536) is 120 MeV heavier 
than the scalar 



TABLE I: The static \ - \ mass splittings AM for n f = 2 
sea quarks [l9|| for different hopping parameters k at f3 = 5.6. 
The numbers in the last column are subject to an additional 
5 % overall scale uncertainty. 





r /a 


mp/mv 


AMr 


AM/MeV 


0.1560 


5.11(3) 


0.834( 3) 


1.16 (9) 


465( 35) 


0.1565 


5.28(5) 


0.813( 9) 


1.15(11) 


460( 45) 


0.1570 


5.48(7) 


0.763( 6) 


1.10(13) 


440 ( 50) 


0.1575 


5.89(3) 


0.704( 5) 


1.24(12) 


495 ( 50) 


0.1580 


6.23(6) 


0.574(13) 


1.08(24) 


430(110) 



Results with nf = 2 mass-degenerate flavours of Wil- 
son sea quarks have been obtained by the SESAM Collab- 
oration [I3 (circles) on slightly finer lattices. The data of 
Table[I]do not exhibit any visible light quark mass depen- 
dence. The interpolated values have been obtained from 
a linear fit in m 2 P (with tiny slope) and the errors of the 
interpolation, that have been conservatively estimated by 
varying fit range and functional form, are dominated by 
systematics. The lattice spacing a obtained from the phe- 
nomenological value a/ro ss ax 400 MeV is in agreement 



with the one obtained from m p a, within errors |19j. The 
experimental lines in the figure correspond to masses, 
relative to the pseudoscalar state, a somewhat arbitrary 
choice since vector and pseudoscalar will be degenerate in 
the static limit. The rationale behind this is that Ref. 
assumes the 1 + - + splitting to be identical to the 1~ 
0~ splitting. The inclusion of sea quarks seems to result 
in the slightly increased value, AM = 468 (43) (24) MeV 
for the s quark system. We do not expect a lattice spac- 
ing dependence of this number in excess of the statistical 
uncertainty, based on the quenched experience. 



IV. FINITE MASS CORRECTIONS 

Effects of the finite charm quark mass have only been 
investigated in the quenched approximation. In particu- 
lar three studies exist: two using lattice NRQCD [20H2H 
to order 1/m 2 and 1/m 3 , respectively (the leading cor- 
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rections are of order a s /m in both cases), and one using 
relativistic charm quarks [2^. Both NRQCD results are 
consistent with each other. The study of Hein et al. [2£j 
has been performed at f3 = 5.7 and j3 = 6.2 for the B s 
and Bd families and at (3 — 5.7 for the D s . The rela- 
tivistic study has been made at f3 = 6.0 and (3 = 6.2. 
In the latter case we refrain from citing values for the 
B meson since the extrapolation of results obtained for 
heavy quark masses much lighter than the b is not fully 
under control. In none of these cases statistically signif- 
icant lattice spacing effects have been observed and we 
display the results for the + - CP splittings in Table ITT1 

TABLE II: The + - 0~ mass splitting in the heavy-light sys- 
tem for two sea quarks in the static limit and in the quenched 
approximation, for the B and D systems in NRQCD [2(1 ] and 
for the D system with relativistic quarks [22J. The errors do 
not include uncertainties in the overall scale which we esti- 
mate to be about 5 % for n/ = 2. All numbers are in units of 
MeV. 





n f = 2 




n f 


= 






static 


static 


NRQCD 


NRQCD 


relativ. 








h = b 


h = c 


h — c 


hs 


468(43) 


384 (50) 


345(55) 


465(50) 


495(25) 


hd 


472(85) 


299(114) 


370(50) 




465(35) 



Note that while the NRQCD results for the B systems 
agree with the respective n/ = static limits, the split- 
ting is enhanced in the D system, in agreement with the 
fully relativistic calculation. The relativistic D s split- 
ting is bigger by as much as (29 ± 16) % with respect 
to the static limit. If we assume a similar increase for 
the case with sea quarks we would expect a splitting of 
600(110) MeV for the D s system, yielding the predic- 
tion m(Dj ) = 2.57(11) GeV. The potential model of 
Ref. [f| predicts 2.48 GeV, while the quenched results 
are 2.44(5) GeV (NRQCD) and 2.47(3) GeV (relativis- 
tic D quark), all significantly bigger than the candidate's 
mass of 2.32 GeV. The quenched lattice results for the D s 
system also indicate a tiny l' + - 1+ splitting, suggesting 
that the 1 + state should be heavier than 2.46 GeV. 



V. SUMMARY 

We calculate a scalar-pseudoscalar splitting of AM = 
468(43) (24) MeV in the static limit for rif = 2 sea quarks, 
significantly larger than the value 338 MeV suggested by 
a heavy quark constituent quark model pl| and larger 
than the quenched QCD value AM = 384(50)(20) MeV. 



We also report a significant finite charm mass correc- 
tion that casts doubt onto naive generalisations to the 
B system. Lattice predictions on the masses are consis- 
tent with the quark model of Ref. and incompatible 
with the new state observed by BaBar and CLEO. We 
conclude that the D^j(2317) might receive a large DK 
component: the physics of this heavy scalar might indeed 
resemble elements of that governing the /n/ao(980) sys- 
tem. If this is the case then the masses of the up and 
down quarks will play a major role and simulations with 
non-mass-degenerate sea quarks are required. 

Unfortunately, on the lattice the possibility of four- 
quark states has so far only been addressed in the 
static limit where attraction was reported in some 
channels |23j | . In view of the new experimental candidate 
quenched simulations of relativistic four quark molecules 
are urgent. To understand the exact nature of the new 
state not only the spectrum but also predictions of decay 
rates are required. While lattice calculations of strong 
decays are unfeasible, a study of electro-magnetic decay 
rates is a possibility. 

Notes added in proof 

The discovery of the two D s mesons has also been con- 
firmed by the Belle Collaboration [24[. 

A new lattice study of D s mesons by the UKQCD Col- 
laboration has appeared recently [25j and a paper by 
Terasaki |26( on the new D s mesons was submitted to 
the preprint server only one day after this article. 

In view of the possibility of similar states in the B s 
spectrum it appears worthwhile to mention that the 
static rif = 2 lattice results presented here imply that 
the scalar quark model B s meson should have a mass 
of 5837(43) (24) MeV, with additional 1/m corrections of 
order 40 MeV, possibly upwards, based on the D s expe- 
rience. This has to be compared with the BK threshold 
of about 5775 MeV. 
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